Electron spin coherence in semiconductors: Considerations for a spin-based solid state 

quantum computer architecture 
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We theoretically consider coherence times for spins in two quantum computer architectures, where 
the qubit is the spin of an electron bound to a P donor impurity in Si or within a GaAs quantum 
dot. We show that low temperature decoherence is dominated by spin-spin interactions, through 
spectral diffusion and dipolar flip-flop mechanisms. These contributions lead to 1 — f 00 /is calculated 
spin coherence times for a wide range of parameters, much higher than former estimates based on 
Tj measurements. 
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Solid state quantum computation based on the intrin- 
sic two level dynamics of electron spin in semiconductors 
has attracted widespread attention because the enormous 
resources of conventional electronics can in principle be 
applied to develop a scalable quantum computer (QC). 
In this context, solid state spins in an applied magnetic 
field are attractive qubit candidates since they comprise 
perfect two-level systems with potentially long coherence 
times. Nuclear spins of phosphorus impurities in sili- 
con are potential qubits because they are well isolated 
from other degrees of freedom [1]. But for such a quan- 
tum computer to work, precise electronic control of single 
nuclear spins must be achieved, a rather daunting task. 
Electron spins are much easier to control. In that re- 
spect, electrons bound to phosphorus impurities in silicon 
[2] and gallium arsenide quantum dots [3] are promising 
qubit candidates. However, for quantum computation to 
be fault-tolerant these electron spins must be coherent 
for at least 10 4 elementary quantum operations [4], which 
imposes a severe constraint since very long spin coherence 
times would be needed. In this letter we address the im- 
portant question of principle involving the fundamental 
upper bound on electron spin coherence time in proposed 
semiconductor QC architectures. Since such a device can 
only work at low temperatures (kgT <C Ez, where Ez 
is the Zeeman splitting of the spins in an applied mag- 
netic field B), we show that the dominant decoherence 
contribution comes from the unavoidable spin-spin inter- 
actions with nuclei and other electrons. This leads to 
an unsurmountable upper bound on spin coherence, but 
our calculations indicate that the fault-tolerance crite- 
rion can still be satisfied for a wide range of parameters 
(e.g. quantum dot Fock-Darwin radius Zo , concentration 
of 29 Si isotopes /, etc) defining the QC architecture. 

An essential property of electrons in Si:P and GaAs 
is that their electron spin resonance (ESR) line is inho- 
mogeneously broadened by the hyperfine interaction with 
nuclear spins [5,6]. This effect leads to a drastic difference 
between the precessing magnetization of an ensemble of 



spins as compared to a single spin or a group of them 
with the same Larmor frequency. The former decays in 
a time scale T 2 * which is dominated by the dephasing ef- 
fect of the inhomogeneous field distribution. The latter 
magnetization decays in a time T 2 which is usually many 
orders of magnitude longer [7] . To extract the single spin 
coherence time T% from an ESR experiment one has to 
perform a ir/2 — tt spin echo sequence [8], that forces the 
spins to refocus eliminating the dephasing effect. This 
leads to a convenient definition for spin coherence: Sim- 
ply the time it takes for a spin echo envelope to decay to 
1/e of its initial value. Since spin echoes usually decay 
in a quite different fashion from the simple exponential 
predicted by the Bloch equations, we denote its decay 
time by Tm using a notation consistent with existing lit- 
erature [9]. Spin echo experiments have been performed 
in dilute Si:P [10,11], but never in GaAs. For the lat- 
ter material, claims of "long coherence times" have been 
based on T 2 * measurements only [12], which do not re- 
veal the ultimate limit on spin coherence set by Tm , that 
can in principle be many orders of magnitude longer than 
T 2 * . Certainly inhomogeneous broadening imposes severe 
tuning constraints on one-qubit gates. But fortunately 
one can still have an universal QC built only with two- 
qubit gates [13]. Our former work shows that the spread 
of Zeeman frequencies will only weakly affect two-qubit 
gates, since the exchange interaction is fairly insensitive 
to inhomogeneous fields [14]. Therefore we emphasize 
that the relevant spin coherence time for QC architec- 
tures is, in fact, Tm (not T 2 ), and our theoretical finding 
of truly large values of Tm Qt M s ) compared with the 
measured T 2 * (~ ns) values is quite significant. 

Due to extreme sensitivity requirements, ESR in GaAs 
is usually measured indirectly, for example from the 
changes in magneto-resistivity in a 2DEG [15] or from 
photoluminescence [16]. In these experiments one probes 
ensemble spin properties of moving electrons or recom- 
bination pairs respectively, meaning that one should be 
careful in extrapolating those results to a single localized 
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spin in a quantum computer environment. Alternative 
methods to ESR include Faraday rotation [12], but again 
to study single spin coherence it would be necessary to 
perform an echo sequence, that is yet to be done using 
optical methods. Measuring single spin coherence time in 
semiconductors is a dauntingly difficult task, and there- 
fore the need for theoretical estimates of Tm becomes 
acutely necessary. Here we present the first realistic cal- 
culation of Tm for a GaAs quantum dot (QD) based QC 
architecture (and also for Si:P QC architecture). Our 
values exceed former coherence estimates based on T 2 * 
by three orders of magnitude, and establish quite defini- 
tively that fault tolerant quantum computation should 
be possible in semiconductor QC at low temperatures. 

First we discuss the spin-flip rate T^ 1 with correspond- 
ing energy transfer to the lattice. Any spin-flip process 
contributes to exponential decay of the spin echo signal. 
But bound electrons at low temperatures quite gener- 
ally have extremely long T\ (3> Tm), since for the spin 
to flip, a virtual transition to an excited orbital state 
must take place, quite different from the conduction elec- 
tron case where any momentum relaxation event may 
flip the spin [17] through spin-orbit coupling. The direct 
phonon emission rate for bound-electrons becomes T^ 1 cx 



[n(E z ) + 1]S 5 , where n{E z 



[exp(E z /k B T) - 1}- 



is the Bose occupation number for the emitted phonon 
wave-vector. Feher and Gere measured Tf »4x 1CP 4 
s" 1 for Si:P at B = 0.32 T and T = 1.25 K [18]. For GaAs 
dots with Fock-Darwin radius lo = 30 nm and B=l T one 
obtains T-f 1 w 200 s _1 [19]. These spin-flip rates when 
compared to other decohcrcnce mechanisms considered 
below give negligible contributions, with the possible ex- 
ception of very large dots (Iq > 100 nm), due to the fact 
that Tj -1 <x Iq, a case we do not consider here since rea- 
sonable semiconductor QC architectures are limited by 
an inter-dot distance of 50 nm, necessary for exchange 
gate operations to work [3]. 

We now describe the spin decohcrcnce mechanisms of 
importance to QC architectures. Two electron spins may 
flip-flop due to their dipolar interaction. But this event 
is limited by inhomogeneous broadening, because only 
spins with the same Zceman splitting can satisfy energy 
conservation in this process [20] . The rate is given by [9] 



1 tl x (Aw 2 ) 
T, * ^ X — 



0.33 



(Aw 2 



(1) 
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Here f(u>o) is the inhomogeneously broadened Gaussian 
line-shape, which plays the role of the density of states 
per unit frequency, (Aw)/ the inhomogeneous line-width, 
while (Aw 2 ) /9 is the flip-flop contribution to the second 
moment [8], which is of the order of the transition ma- 
trix element squared. It is important to mention that 
the qubit-qubit dipolar interaction can in principle be 
included in the Hamiltonian responsible for the quantum 
algorithm [21], and also may be eliminated using mag- 
netic resonance techniques [8]. But here we have two 



reasons to include it in our calculation: First we do not 
expect forthcoming Tm measurements to be free from 
these effects. Second, since correcting dipolar coupling 
implies additional overhead in QC design, it is interest- 
ing to access the amount of error involved by ignoring its 
presence. Another spin-spin mechanism is spectral diffu- 
sion, that happens when the electron spins that generate 
the echo pulse are subject to fluctuating dipolar or hy- 
perfine fields generated by nuclear spins. A stochastic 
theory of this effect can be formulated by treating the 
electron spin Larmor frequency w as a random variable, 
and calculating the echo envelope amplitude M(2r) using 
an ensemble average [22]. The result is 



M(2r) w M(0)exp [-{RS\/ 12) (2r) 



(2) 



where r is the time interval between the it/ 2 and ir pulses, 
and R the local field relaxation to a Gaussian probability 
distribution with width 8 a- Hence the spectral diffusion 
(SD) rate is given by T<7^ = (RS 2 J 12) 1 / 3 . The deco- 
herence rate follows from all the above contributions by 
solving the cubic equation for T M X : (T^ 1 + TJ x )Tm + 
( t sd t m) 3 = 1, and often one finds it to be dominated 
by a single mechanism. (We give our calculated values 
for Tm in Si:P and GaAs quantum dot system in Table 
1, where the dominant decoherence mechanism in each 
case is also listed.) 

Now we estimate Tm for Si:P. For the dipolar mech- 
anism, assuming a qubit separation d = 10 nm we get 
Tj 1 = 1 x 10 3 s" 1 and 2 x 10 3 s" 1 when the qubits 
are arranged in ID and 2D lattices respectively. This is 
the case for 29 Si natural abundance (/ = 4.67%) that 
leads to an experimental line- width of 2.5 G [5]. How- 
ever, for isotopically pure 28 Si (which has zero spin), we 
expect no inhomogeneous broadening [(Aw)/ ~ (Aw 2 ) 
in Eq. (1)]. In that case, TJ 1 will be much higher: 
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(2(Aw 2 )/tt) 1 / 2 = 3 x 10 5 s-\ 4 x 10 5 



domi- 



nated by qubit dipolar coupling. The SD rate was mea- 
sured for natural silicon by Chiba and Hirai [11]. By as- 
suming that fluctuating dipolar fields of 29 Si nuclear spins 
caused Gaussian spectral diffusion they formulated a the- 
ory to calculate the coefficient R8\ in (2), which agreed 
within order of magnitude with the measured Tm (note 
that since their P concentration was low TJ 1 was negli- 
gible). Therefore we assume their experimental value of 
Tg D — 3.3 x 10 3 s _1 as a reliable estimate for the SD 
rate for a single P spin in a Si matrix, and by adding the 
dipolar flip-flop rate we are able to get a phcnomcnolog- 
ical estimate for the decoherence rate when the P spins 
are arranged in a quantum computer geometry, as op- 
posed to a dilute randomly doped sample. Hence we get 
T M l «4x 10 3 s _1 in both ID and 2D geometries, noting 
that here both spin-spin mechanisms contribute with the 
same order of magnitude. An interesting consequence of 
the interplay of these two mechanisms is that T M X dis- 
plays a minimum as a function of the 29 Si fraction /, 
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an effect not yet noted in the literature (Fig. 1). The 
Tgp contribution is proportional to J 2 / 3 , since the prob- 
ability of finding a pair of 29 Si [analogous to pij in (10) 
below] is proportional to / 2 , and then we take a cubic 
root due to the non-exponential decay in (2). However, 
the dipolar flip-flop rate is decreased when we increase 
/: TJ 1 oc (Auj)J 1 oc / _1 / 2 . This happens because the 
29 Si is the source of hyperfine broadening and the mean 
square deviation of the hyperfine field is proportional to 
/. It would be interesting to fit these results to measure- 
ments on isotopically purified samples; to date there is 
only one measurement, with / = 0.12%, in which T^ 1 
decreased by a factor of 2 [10] (the P concentration was 
the same as in Fig. 1, c = 4 x 10 16 cm -3 ). This is con- 
sistent with our result. For the case of a 2D Si:P QC 
we find that T^ 1 will attain a minimum when / w 2%, 
suggesting that natural silicon is a good choice for QC 
architectures. 

Turning to the GaAs-QD, we note that inhomogeneous 
broadening should be much stronger: We estimate a line- 
width of at least 50 G due to hyperfine interactions, since 
all nuclei have spin 3/2 (S impurities in GaAs have 500 G 
of broadening [6]). This together with the fact that the 
qubits will be much further apart (d = 50 nm) leads to 
TJ 1 < 10~ 3 s" 1 (Tf 1 oc d- 6 ). We now calculate the SD 
rate for these spins. NMR in GaAs reveals a composition 
of 75 As (50%), 69 Ga (30.2%), and n Ga (19.8%) [23]. The 
homogeneous line- width of the 75 As system (fee lattice 
with a = 5.65 A) is 5 B = 5.56 x 10 3 s _1 [24], indicating 
that these nuclei are flip-flopping every 300 us. Since 
the hyperfine interaction depends on the position of each 
nucleus, the electron feels a different field if a pair of 
nuclei is up-down as opposed to down-up (here we will 
neglect the contribution of nuclei outside the QD wave 
function; they would produce a dipolar field which shifts 
the electron frequency by at most 5u) <~ 7e7nV z o ~ 20 



< S r 



t; 



200 s , the intrinsic line- width due 



to QD spin relaxation mentioned above; this is certainly 
not the case in Si:P, where the nuclei outside the donor 
wave function give the dominant contribution [11]). To 
calculate the spectral diffusion rate we consider a three 
spin Hamiltonian H = H a + V , where 



Hq — hoJiIi z S z + huJjIj z S z , 

V = nbij(I i+ Ij- + h-Ij+) - 4:hbijI iz I jz , 
1 , . 1 — 3 cos 2 8^ 



bij 



Lo t = Lu h exp 



R 3 

Xf + Y( 
I 2 



COS 7T 



(3) 
(4) 

(5) 
(6) 



Here H$ is the diagonal part of the hyperfine interac- 
tion (the off-diagonal component can be neglected since 
electron-nucleus flip-flop is forbidden by energy conser- 
vation), while V is the dipolar coupling between two nu- 
clei. Ii is the spin operator for the nucleus located at 



R ; = (Xi,Yi,Zi), 7n = 4.58 x 10 3 (s G)" 1 the gyro- 
magnetic ratio for 75 As, Rij the distance between the 
two flip-flopping nuclei and % the angle between this 
vector and the magnetic field. The hyperfine frequency 
is given by 

8tt llT ,„,,2 16 7 e 7 n ?ian , „ 
= y7e7n?i l*(0)| 2 = y ^ °- d(As), (7) 

where d(As) = |u'(0)| 2 is the electronic density on the 
nuclei and lo the Fock-Darwin length. The electron fre- 
quency change due to a nuclear flip-flop is then Qij = 
uji — u)j, while the energy changes by KQij/2 in this pro- 
cess. Since \bij\ < |Qjj| we can apply perturbation theory, 
and the rate becomes [11] 



r« = 27r6 2 j5B (r! y /2), 



2¥5 B eXP V 2<5| 



(8) 
(9) 



The line-width 8b of the nuclear system enters to guar- 
antee energy conservation: the change in the Zeeman 
energy of the electron is compensated by the spin-spin 
interaction of the nuclear system. The SD rate can now 
be estimated by summing over all pairs ij, as long as 
they can flip-flop: 
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!a{oj) = 



■ exp 



Pij(T), (10) 
(11) 



/2ir6 r * " V 2^/' 
Assuming the nuclei are in thermal equilibrium we get 

cosh3x 



Pij(T) 



1 + 



2 cosh 2x + 3 cosh x + 2 



(12) 



with x = h~f n B/k B T. For B = 1 T and T > 1 mK 
we assume the high T limit w 7/8. By performing 
the sum (10) numerically, we estimate Tgf, ~ 10 4 s _1 , 
that shows for QD's with Z ^ 100 nm spectral diffusion 
due to hyperfine field fluctuation dominates. This rate 
depends on B field intensity (Fig. 2), since lo decreases 
when B is increased, and also varies by about a factor of 
2 when the tilting angle with z direction is changed [see 
Eq. (5)]. 

Finally, we wish to comment on the validity of the ap- 
proximations employed here. The SD decay given by Eq. 
(2) is only valid for Rt < 1 [22], and also Eq. (10) as- 
sumes a generalization of this two parameter model to 
several parameters r^, flij. Certainly a more rigorous 
theory of SD due to nuclear spins needs to be developed 
if one wants to go beyond the order of magnitude esti- 
mates given here. In particular, a rigorous theory for the 
nuclear flip-flop rate [Eq. (8)] is mandatory for a precise 
description of this phenomena. Our estimated coherence 
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times, summarized in Table 1, should be compared with 
the longest gating time in the corresponding QC archi- 
tecture. For the parameters chosen here, the exchange 
time will be tj ~ h/O.l meV ~ 1 ps, with a typical 
exchange coupling [3] of 0.1 meV. A single qubit rota- 
tion (Rabi flop) can be done with an ESR field 20 times 
smaller than the applied field, that leads to a ir/2 rota- 
tion time t r ~ 2O/7B - O.f ns, 0.5 ns at B = 1 T for 
Si:P, GaAs-QD respectively. These time scales lead to a 
quality factor Tm/tr > 10 4 for both architectures. Also, 
Tm/tj ~ f0 6 - f0 8 (since tj ~ 1 ps and T M ~ 1 - 100 
/its), which implies a large number 10 4 ) of coherent 
gate operations allowing convenient fault tolerant com- 
putation well within the currently estimated (10 ) error 
correction scheme. Hence the electron spin in semicon- 
ductors is confirmed as a competitive qubit candidate, 
with the effective low temperature upper bound on the 
coherence time given by 1 — 100 jus under quite general 
conditions. The authors acknowledge discussions with J. 
Fabian, X. Hu, A. Kaminski, B. Koillcr, and I. Zutic. 
This work is supported by ARDA, LPS, US-ONR, and 
NSF. 
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TABLE I. Coherence times Tm, quality factors, and dom- 
inant decoherence mechanisms for three quantum computer 
architectures. We assume the qubits are disposed in a 2D 
square lattice of side 10 nm and 50 nm for the case of Si:P 
and GaAs-QD respectively. For the GaAs-QD architecture, 
Fock-Darwin radius is assumed 30 nm. 



4 



10 



03 
O 

" 1 

CD 

cd 
i — 

o 
o 

S 0.1 h-'.' 



CD 

o 
o 

CD 

Q0.01 



0.005 



2D Square latt. 
• Random doping 



4X10 16 cm" 3 



j aI6 -3 

10 cm . 



Natural Si 



0.01 



0.1 



Amount of Si [%] 



FIG. 1. Depicts the competition between the dipolar 
flip-flop rate and spectral diffusion for a bound-electron spin 
in Si:P, leading to a minimum in the decoherence rate as 
a function of the 29 Si fraction. We show calculations for a 
2D quantum computer architecture with qubit separation of 
lOnm, and for the most common experimental situation of 
random phosphorus doping, with concentrations 1 — 4 x 10 16 
cm" 3 . 
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FIG. 2. Calculated spectral diffusion rate as a function of 
magnetic field and quantum dot Fock-Darwin radius lo (in- 
set). For the B field plot the dot transverse confinement 
length is set to 50 nm. The spectral diffusion mechanism com- 
pletely dominates the decoherence rate for small dots (lo <100 
nm). 
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